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Abstract 

Using the CLEO detector at the Cornell Electron Storage Ring we have studied the angu- 
lar distributions in the decay — > Ke + u e . By performing a four-dimensional maximum like- 
lihood fit, we extract the form factor ratio, R = /2//1, and the pole mass, which determines 
the shape of the form factors, M po i e . They are found to be —0.31 ± 0.05 s t a t ± 0.04 SJ/S and 
(2.13 ± 0.07 s tat ± O.lOsys) GeV/c 2 respectively. These results correspond to the following value 
of the decay asymmetry parameter a\ c = —0.85 ± 0.03 s t a t ± 0.02 sys , for (q 2 ) = 0.67 (GeV/c 2 ) 2 . 
We search for CP violation in the angular distributions of the decay and find no evidence for CP 
violation: Aa c = {a c _ a ^ c = 0.01 ± 0m sta t ± 0M sys ± 0.02^ A , where the third error is from the 

uncertainty in the world average of the CP violating parameter, Aa, for A — > p7r~. All results 
presented in this paper are preliminary. 



* Submitted to the International Europhysics Conference on High Energy Physics, July 2003, Aachen 
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Charm semileptonic decays allow a measurement of the form factors which parameterize 
the hadronic current because the Cabbibo-Kobayashi-Masawa (CKM) matrix element \V CS \ is 
known from unitarity pj. Within heavy quark effective theory (HQET) 0, A-type baryons 
are more straightforward to treat than mesons as they consist of a heavy quark and a spin 
zero light diquark. This simplicity allows for more reliable predictions concerning heavy 
quark to light quark transitions |3j, |4| than is the case for mesons. A measurement of the 
form factors in A+ — > Ae + u e will help the future determination of the matrix element \V u b\ 
using A° decays since HQET relates the form factors in A+ decay to those governing A° 
semileptonic decays. 

In the limit of negligible lepton mass, the semileptonic decay of a charmed baryon 
(l/2 + — > l/2 + ) is usually parameterized in terms of four form factors: two axial form factors 
F A and F^ and two vector form factors F, and _F 2 V ". These form factors are functions of q 2 , 
the invariant mass squared of the virtual W + . In the zero lepton mass approximation, the 
decay may be described in terms of helicity amplitudes H Xa x w = H\ A \ W + H\ A \ w -> where Aa 
and \y/ are the helicities of the A and W + . The helicity amplitudes are related to the form 
factors by 0| 



where Q± 

the parity relations H_ 



V 2 u 


= jQl[(M Kc +M A )FY-q 2 FV], 
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hY}^L. In terms of the helicity amplitudes, the decay 



angular distribution can be written as 







dT - B(\ -) 1 G<2p W I 2 q2p 

~ dq 2 d cos 8 A d cos 8 w d X ~ ( ^ W h~(2^y\ cs| 24MX X 

3 3 
{-(1 - cos^) 2 |^ii! 2 (l + a A cos8 A ) + zi l + cos 0w) 2 |^T_i_ 1 | 2 (l - a A cos# A ) ( 2 ) 

o 2 o 2 

3 

+ - sin 2 w [|.ffi | 2 (l + a A cos# A ) + \H_^ \ 2 {1 - a A cos# A )] 
3 

H i=a A cos x sin 0^/ sm^A^l — cos9 w )Re(H_i Hl 1 ) + (1 + cos 9w)Re(Hi H^_ i_ 1 )]}, 

2y2 2 2 2 2 

where Gf is the Fermi coupling constant, V cs is the CKM matrix element, P is the A 
momentum in the A+ rest frame, 8 A is the angle between the momentum vector of the 
proton in the A rest frame and the A momentum in the A+ rest frame, 8\y is the angle 
between the momentum vector of the positron in the W + rest frame and the A momentum 
in the A+ rest frame, x is the angle between the decay planes of the A and W + , and a A is 
the A — > p-K~ decay asymmetry parameter measured to be 0.642 ± 0.013 yj. 

Within the framework of HQET the heavy flavor and spin symmetries imply relations 
among the form factors which reduce their number to one when the decay involves only heavy 
quarks. In this Letter we follow Ref. 0] where it is argued that the HQET formalism may 
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also work for A+ — > Ae + u e . Treating the c quark as a heavy quark, two independent form 
factors fi and f2 are required to describe the hadronic current. The relationships between 
these form factors and the standard form factors are F^(q 2 ) = —F^(q 2 ) = fi(q 2 ) + jf^-f2(q 2 ) 

and F^iq 2 ) = —F^(q 2 ) = j^t^i.^ 2 )- In general / 2 is expected to be negative and smaller 
in magnitude than f\. If the strange quark is treated as heavy, f% is zero. 

In order to extract the form factor ratio R = /2//1 from a fit to T an assumption must 
be made about the q 2 dependence of the form factors. We follow the model of Korner and 
Kramer (KK) 0] who use the dipole form f(q 2 ) = ^J^J^ )2 (1 - |p^) 2 , where the pole 

mass is chosen to be m^.^-) = 2.11 GeV/c 2 . We will perform a simultaneous fit for the 
form factor ratio and the pole mass. 

The data sample used in this study was collected with the CLEO II [6] and the up- 
graded CLEO II. V [?J detector operating at the Cornell Electron Storage Ring (CESR). 
The integrated luminosity consists of 13.7 fb -1 taken at and just below the T(4S) reso- 
nance, corresponding to approximately 18 x 10 6 e + e~ — *> cc events. Throughout this paper 
charge conjugate states are implicitely included, unless otherwise is indicated, and we use 
the symbol e to denote an electron or positron. 

The analysis technique is an extension of previous work |j| . We search for the decay 
A+ — > Ae + v e in e + e~ — > cc events by detecting a Ae + pair with invariant mass in range 
m\ e + < m\ c |8||. All tracks are required to come from the region of the event vertex. To 
reduce the background from B decays, we require R2 = H 2 /H Q > 0.2, where Hi are Fox- 
Wolfram event shape variables llO]. Positrons are identified using a likelihood function which 
incorporates information from the calorimeter and dE/dx systems. The minimum allowed 
momentum for positron candidates is 0.7 GeV/c, as the positron fake rates are much higher 
in the lower momentum range. Positrons are required to have been detected in the region 
I cos#| < 0.7, where 9 is the angle between the positron momentum and the beam line. 

The A is reconstructed through its decay to pir~ . We require the point of intersection of 
the two charged tracks, measured in the r — <p plane, to be greater than 5 mm away from 
the primary vertex. In addition, we require the sum of the p and 7r~ momentum vectors to 
extrapolate back to the beam line. The dE/dx measurement of the proton is required to 
be consistent with the expected value. We reject combinations which satisfy interpretation 
as a Kg. Finally, we require the momentum of the pn~ pair to be greater than 0.8 GeV/c 
in order to reduce combinatoric background. These A candidates are then combined with 
Right Sign (RS) tracks consistent with positrons, and the sum of the A and e + momenta is 
required to be greater than 1.4 GeV/c, in order to reduce the background from B decays. 

The number of events passing these requirements is 4060, of which 123 ± 12 are consistent 
with fake A background, 338 ± 67 with S c — > Se + z/ feedthrough and 398 ± 58 are consistent 
with the e fake background. The sidebands of the p7r~ invariant mass distribution are used 
to estimate the fake A. The background from S c — > Se + u decays is estimated using the 



result of a previous CLEO analysis |llj |. The e fake background is estimated using the 
Wrong Sign (WS) h + A (no charge conjugation is implied here) data sample, where h + A are 
WS combinations passing the analysis selection criteria. The h + tracks in this sample are 
mostly fakes as there are few processes contributing e + A pairs after the selection criteria are 
applied. We expect that the estimate of the normalization and the momentum spectrum 
of the e fake background from the h + A sample is accurate for the following reasons. The 
probability to form a h + A or h~A pair is approximately equal because the net charge of the 
event is zero. The baryon number and the (sud) quark content of the A imply (1) that, due 
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to baryon conservation, a A is more likely to be produced with an antiproton in WS rather 
than RS combinations and (2) that if a A is produced as a result of a ss quark pair creation, 
such events have a higher fraction of kaons in RS rather than WS combinations. In addition 
antiprotons and kaons have higher e fake rates. By using only one charge conjugate state 
(h + A) and by excluding the momentum region where the e fake rate from kaons is high by 
requiring \p e \ > 0.7 GeV/c 2 , differences between the momentum spectra and particle species 
of hadronic tracks between h + A and h + A are minimised. Differences that remain are second 
order and are accounted for as a systematic uncertainty on the final result. 

Calculating kinematic variables requires knowledge of the A+ momentum which is un- 
known due to the undetected neutrino. The direction of the A+ is approximated based on 
the information provided by the thrust axis of the event and the kinematic constraints of 
the decay. The magnitude of the A+ momentum is obtained as a weighted average of the 
roots of the quadratic equation p A 2 c = (p\ + Pe + Vv) 2 ■ The weights are assigned based on 
the fragmentation function of A+. After the A+ momentum is estimated, the four kinematic 
variables are obtained by working in the A+ center-of-mass frame. 

Using t = q 2 /q 2 nax , cos^a, cos 0w, and we perform a four-dimensional maximum likeli- 
hood fit in a manner similar to Ref. 12] . The technique enables a multidimensional likelihood 
fit to be performed to variables modified by experimental acceptance and resolution and is 
necessary for this analysis due to the substantial smearing of the kinematic variables. The 
essence of the method is to determine the probability density function by using the popula- 
tion of appropriately weighted Monte Carlo events in the four-dimensional kinematic space. 
This is accomplished by generating one high statistics sample of Monte Carlo events with 
known values of the form factor ratio R, the pole mass M po i e and corresponding known 
values of the four kinematic variables t, cos9\, cos 9w, and \ f° r eacn event. The generated 
events are then processed through the full detector simulation, off-line analysis programs, 
and selection criteria. Using the generated kinematic variables, the accepted Monte Carlo 
events are weighted by the ratio of the decay distribution for the trial values of R and M po i e 
to that of the generated distribution. The accepted Monte Carlo events are now, therefore, 
distributed according to the probability density corresponding to the trial values of R and 
Mp i e . By such weighting, a likelihood may be evaluated for each data event for different 
values of the form factor ratio and the pole mass, and a fit performed. 

The fit is unbinned in cos#a, cos 6^ and x an d takes into account the correlations among 
these variables. The probability for each event is determined by sampling this distribution 
using a search volume around each data point. The volume size is chosen so that the 
systematic effect from finite search volumes is small and the required number of Monte 
Carlo events is not prohibitively high. We bin the fit in t for the following reason. In the 
course of this analysis we observed an inconsistency between datasets for one subsample of 
data for 0.8 < t < 1.0 (a total of 237 evnets). The projections in the other three variables 
for this subsample are consistent with the remainder of the data. To avoid bias we wish to 
include these events in the fit for three variables, but not for t. This is not possible if the 
fit is unbinned in 4 dimensions. Instead we include the events in the 3 dimensional fit but, 
by virtue of perfoming a binned fit in t, we exclude their contribution to the calculation of 
the likelihood in t. A systematic uncertainty is assigned by varying the excluded region. 

The A+ — > Ae + v sample has signal to background in the approximate ratio 4:1. Back- 
ground is incorporated into the fitting technique by constructing the log-likelihood function 
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FIG. 1: Projections of the data (points with error bars) and the fit (solid histogram) onto t, cos 9\, 
cos6>v^ and x- The dashed lines show the sum of the background distributions. 
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The first term is the sum over the number of events, where N events is the number of events 
in the signal region, P s and Pj B are the probabilities for the ith event to be signal and 
belong to the jth background component, respectively, Y l s is the signal shape modeled from 
signal Monte Carlo, and H B is the background shape modeled by a sample of events for 
the jth background component. The second term is the sum over bins of t, where n l observed 
is the number of events in the zth bin, n l estimated is our prediction of that number, which is 
dependent on the trial values for R and M po i e and is estimated similarly to T s and r*- B . 

The e fake background is modeled by the fake positron data sample. Feedthrough back- 
ground from S c — > Se + z/ is modeled by the Monte Carlo sample which is generated according 
to the HQET-consistent KK model. Fake A background is modeled using the data events 
in the sidebands of the pir~ invariant mass distribution. 

Using the above method we perform a simultaneous fit for the form factor ratio and 
the pole mass and find R = -0.31 ± 0.05 and M pole = (2.13 ± 0.07) GeV/c 2 , where the 
uncertainties are statistical. This is our main result. Figures ^ and El show the t, cos 6a, 
cos 9w and x projections for the data and for the fit. 

We have considered the following sources of systematic uncertainty and give our estimate 
of their magnitude in parentheses for R and M po i e . The uncertainty associated with the 
size of the search volume is measured from a statistical experiment in which a set of mock 
data samples, including signal and all background components, was fit in the same way as 
the data (0.006, 0.048). The uncertainty due to the limited size of the Monte Carlo sample 
is estimated by dividing the Monte Carlo sample into four independent equal samples and 
repeating the fit (0.007, 0.012). The uncertainty due to background normalizations is de- 
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cos8 A c°s6 w X 

FIG. 2: Projections of the data (points with error bars) and the fit (solid histogram) onto cos8\, 
cosOw and x f° r two t regions. The plots labeled (a), (b) and (c) are for t < 0.5; (d), (e) and (f) 
are for t > 0.5. The dashed lines show the sum of the background distributions. 



termined by varying estimated number of the background events in the signal region by one 
standard deviation for each type of background separately (0.023, 0.024). The uncertainty 
associated with the modeling of the background shapes, including uncertainties originating 
from the unknown form factor ratio and pole mass for the decay S c — > Se + z/, are estimated 
by varying these shapes and by using different background samples (0.024, 0.049). The 
uncertainty from random e + A pairs from in qq and T(AS) — ► BB events is estimated by re- 
peating the fit with and without correcting for this background contributions (0.013, 0.038). 
The modes A+ — > AXe + z/, where X represents unobserved decay products, have never been 
observed. The current upper limit is i3(A+ — ► AXe + z/)/£>(A+ — > Ae + u) < 0.15, where 
X 7^ 0, at 90% confidence level The uncertainty due to the possible presence of these 
modes is estimated from a series of fits, each with an additional background component rep- 
resenting possible A+ — ► AXe + v modes. The normalizations of the additional components 
are allowed to float in the fits. By performing fits with a variety of A+ — > AXe + v modes, and 
combining the results, we obtain an estimate of the uncertainty in the final result due to the 
possible presence of A+ — > AXe + u (0.020, 0.060). The uncertainty associated with the A+ 
fragmentation function is estimated by varying this function (0.003, 0.002). The uncertainty 
associated with Monte Carlo modeling of slow pions from A decay is obtained by varying 
this efficiency according to our understanding of the CLEO detector (0.004, 0.003). Finally, 
we account for the effect of excluding 0.8 < t < 1.0 by varying the size of the excluded 
region (0.000; 0.014). Adding all sources of systematic uncertainty in quadrature, our final 
result is R = -0.31±0.05 siat ±0.04 sys and M pole = (2.13±0.07 stat ±0.10 S j /s ) GeV/c 2 . We also 
find R = -0.32±0.04 stat ±0.03 SJ/s for a fit for the form factor ratio with M pole = 2.11 GeV/c 2 
fixed. 

Using the value of R and M po i e obtained in the simultaneous fit and the HQET-consistent 
KK model, the mean value of the decay asymmetry parameter of A^ — > Ae + u e averaged 
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over the charge conjugate states is calculated to be oa c = —0.85 ± 0.03 sta t ± 0.02 s?/s , for 
( g 2) = o.67 (GeV/c 2 ) 2 . 

In the Standard Model CP violation is expected to be very small for semileptonic decays. 
In the KK model, the shape of the decay rate distribution in the 4D space of t, cos8\, 
cos9w and x is governed by R and M po [ e only. We determine ot\ c a^ and ckj by repeating 
the simultaneous fit to R and M po [ e for each charge conjugate state separately. We find 
«a c «a = — 0.561 ±0.026 siai and a j — ~~ 0.535 ±0.024 siai . Following jl3j] and by extension 

we define the CP violating asymmetry of the A+ as A\ c = ( QA C _Jl c ) • From the measurement 
of the CP asymmetry in the product of ax c <y\, and using the relation 

A Ac +A A + 0(A 2 Ac ,A 2 A ), (4) 



we obtain A\ c = 0.01 ± 0.03 sta t ± 0.01 sys ± 0.02_4 A , where in the systematic uncertainty we 
have taken into account the correlations among the systematic uncertainties for the charge 
conjugate states and the third error is from the uncertainty in the world average of the CP 
violating parameter, Aa, for A — > pir~ [l|. 

In conclusion, using a four-dimensional maximum likelihood fit the angular distributions 
of A+ — > Ae + v e have been studied and the form factor ratio R = /2//1 and M po i e are 
found to be -0.31 ± 0.05 stat ± 0.04 sys and (2.13 ± 0.07 stat ± 0.10^,) GeV/c 2 respectively. 
These results correspond to the following value of the decay asymmetry parameter «a c = 
-0.85 ± 0.03 stat ± 0.02 sys , for (q 2 ) = 0.67 (GeV/c 2 ) 2 . We have searched for CP violation 
in the angular distributions of the decay and find no evidence for CP violation: *4a c = 
0.01 ± 0.03 s t a t d= 0.01 S y S ± 0.02a a , where the third error is from the uncertainty in the world 
average of the CP violating parameter, A\, for A — > pn~ . All results presented in this paper 
are preliminary. 
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